Reply to Comment LCK1019 by Prokof'ev and Svistunov 
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The purpose of Ref. was to study the geometry of 
real vortex loops (VLs) in 3D superconductors (3DSC) 
and superfluids (3DSF), and its connection to critical- 
ity. 3DSC are dual to 3DSF, and vice versa [J]. VLs 
originate with transverse fluctuations in the phase 9 of 
the complex pairing field. A fundamental topological 
constraint is that at every point in a 3DSC or 3DSF, 
one must have V • n = 0; n is the vorticity defined by 
V x (V0) = 2nn. Vortices in 3DSC/3DSF cannot start 
or end inside the system. 

We related the Hausdorff dimension Dh of the VLs to 
the anomalous scaling dimension r\ of corresponding field 
theory. In the case of VL, this field theory is the dual of 
the original theory. The relation is r\ + Dh = 2, where rj 
is defined via the correlation function of the field theory 
for the VL (i.e. the dual), G(x,y) = (0*(x)0(y)}. At 
criticality, it takes the form G(x, y) = f/|x — y\ d ~ 2+r i 
in d dimensions and yields the probability of connect- 
ing two points x and y with any directed continuous 
vortex path 0. It is implicit y|/ that any path en- 
tering in computing G(x, y) must be part of a (closed) 
VL. Thus, G(x,y) = N L (x,y)/N L where V L (x,y) is 
the number of VLs passing through both x and y, 
while Nl is the total number of VLs. A useful quan- 
tity is the corresponding probability P„ (x, y) of connect- 
ing x and y with a connected vortex path of length 
n, which must also be part of a closed VL. Hence, we 
have P„(x,y) = N L (x,y;n)/N L , where N L (x,y;n) is 
the number of VLs in the system passing through both 
x and y, when these two points are separated by paths 
of length n. Clearly, G(x, y) = ^2 n P n (x,y)- We used 
the scaling form P„(x,y) = F(|x - y\/n 1 / DH )/n d / DH 
0, and assumed F(0) + 0, in order to have a finite 
probability of returning to the starting point by going 
around the entire closed VL. (Had we for instance as- 
sumed a form F(t) ~ t < 1, & > 0, we would have 
reached the erroneous conclusion that the probability of 
returning to the starting point by following a strictly 
closed VL path would be zero.) Moreover, the distribu- 
tion function D(n) of closed loops of length n is given 
by D(n) = (l/V) ]T x P„(x,x). At the critical point, 
D(n) ~ n~ a , with a = d/D H + 10. 
The graphs of a high-temperature (HT) expansion of the 
partition function Z of the 3DXY model are closed paths 
that can be identified mathematically with sterically in- 
teracting closed loops 0, 0, The graphs in the HT ex- 
pansion of the two-point correlation function G are open- 
ended paths (OEP) 3]. Thus, the geometry of graphs 
contributing to Z differs from that of the graphs con- 
tributing to G • The HT-graphs contributing to Z and 
G can in any case not be physically identified with VLs. 


Ref. computes Dh of OEPs originating with a HT 
graph expansion of (</>*(x)0(y)) for a \<j>\ -theory |5j,|6J, 
which is in the same universality class as the 3DXY 
model. This entails studying an analog of the above cor- 
relation function, namely |3j,|5j G(x, y) = J2 n P«( x j y) = 
J2 n Zn(x,y)/zn, where z n = S y 2 n (x,y) is the number 
of open-ended paths starting in x and ending up in any 
endpoint y, while z n (x, y) is the number of graphs of 
length n starting in x and ending in y. Importantly, 
z n clearly grows with the number of steps n. We write 
Zj. ~ n°l DH , and fUx,y) = F(|x - y\/n 1 / Dll )/ri d / DH 
0,0- In Refs. one effectively defines a loop dis- 

tribution function D(n) ~ n~ a even for these OEPs, by 
setting x — y = a in P n ; a is a lattice constant. 
Again a — d/Dn + 1 but the factor z n in P ra (x, y) 
now requires, for consistency, that F(t) ~ t . Includ- 
ing the scaling relation proposed in [lj is modified as 
follows: rj + Dh = 2 — where i? determines the asymp- 
totic number z n of open-ended paths of length n at the 
critical point in a system of OEPs yj/. Computing Dh 
of such OEPs, Ref. 5] reports D H = 1.7655(20). Using 
i] = 0.038(4) (appropriate for the field theory of closed 
steric loops) 0], they infer d = 0.1925(29) 0. 
While the definition of closed paths used in Ref. |5j is 
appropriate for OEPs, it is not for VLs. OEPs are clearly 
less compact than VL's, whence their Dh is smaller. 
Moreover, OEPs are unphysical objects in the context 
of VLs in 3DSC, as is correctly implied on pg. 3 of Ref. 
0. Hence, computing Dh for OEPs and relating it to n 
of the dual field theory for VLs in a 3DSC (j,|5j, provides 
no test of the relation r) + D H = 2. In 3DSC and 3DSF, 
$ is fixed to zero by the topological constraint V • n = 0. 
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